home *** CD-ROM | disk | FTP | other *** search
/ The Atari Compendium / The Atari Compendium (Toad Computers) (1994).iso / files / umich / falcon / programm.ing / nt_dsp1.lzh / NT_DSP1.MSA / FLOAT / FPMAC.ASM < prev    next >
Assembly Source File  |  1989-01-24  |  3KB  |  88 lines
  1. ;
  2. ; This program originally available on the Motorola DSP bulletin board.
  3. ; It is provided under a DISCLAIMER OF WARRANTY available from
  4. ; Motorola DSP Operation, 6501 Wm. Cannon Drive W., Austin, Tx., 78735.
  5. ; 
  6. ; Last Update 30 Oct 87   Version 2.1
  7. ;
  8. fpmac   ident   2,1
  9. ;
  10. ; MOTOROLA DSP56000/1 FPLIB - VERSION 2
  11. ;
  12. ; FPMAC - FLOATING POINT MULTIPLY/ACCUMULATION SUBROUTINE
  13. ;
  14. ; Entry points: fmac_xya        R = A + X * Y
  15. ;               fmac_mxya       R = A - X * Y
  16. ;
  17. ;       m = 24 bit mantissa (two's complement, normalized fraction)
  18. ;
  19. ;       e = 14 bit exponent (unsigned integer, biased by +8191)
  20. ;
  21. ; Input variables:
  22. ;
  23. ;   X   x1 = mx  (normalized)
  24. ;       x0 = ex
  25. ;
  26. ;   Y   y1 = my  (normalized)
  27. ;       y0 = ey
  28. ;
  29. ;   A   a2 = sign extension of ma
  30. ;       a1 = ma  (normalized)
  31. ;       a0 = zero
  32. ;
  33. ;       b2 = sign extension of ea (always zero)
  34. ;       b1 = ea
  35. ;       b0 = zero
  36. ;
  37. ; Output variables:
  38. ;
  39. ;   R   a2 = sign extension of mr
  40. ;       a1 = mr  (normalized)
  41. ;       a0 = zero
  42. ;
  43. ;       b2 = sign extension of er (always zero)
  44. ;       b1 = er
  45. ;       b0 = zero
  46. ;
  47. ; Error conditions:     Set CCR L=1 if floating point overflow.  Result
  48. ;                       is set to the maximum floating point value of the
  49. ;                       correct sign.  The CCR L bit remains set until
  50. ;                       cleared by the user.
  51. ;
  52. ;                       Set CCR L=1 if floating point underflow.  Result
  53. ;                       is set to floating point zero.  The CCR L bit
  54. ;                       remains set until cleared by the user.
  55. ;
  56. ; Assumes n0, m0, shift constant table and scaling modes
  57. ; initialized by previous call to the subroutine "fpinit".
  58. ;
  59. ; Alters Data ALU Registers
  60. ;       a2      a1      a0      a
  61. ;       b2      b1      b0      b
  62. ;       x1      x0      y1      y0
  63. ;
  64. ; Alters Address Registers
  65. ;       r0
  66. ;
  67. ; Alters Program Control Registers
  68. ;       pc      sr
  69. ;
  70. ; Uses 1 location on System Stack
  71. ;
  72. fmac_xya move   a,fp_space:fp_temp+1    ;save ma
  73.         move    b,fp_space:fp_temp+2    ;save ea
  74.         jsr     fmpy_xy                 ;multiply x*y first
  75.         move    fp_space:fp_temp+1,x1   ;restore ma
  76.         move    fp_space:fp_temp+2,x0   ;restore ea
  77.         jmp     fadd_xa                 ;add product to a
  78. ;
  79. fmac_mxya move  a,fp_space:fp_temp+1    ;save ma
  80.         move    b,fp_space:fp_temp+2    ;save ea
  81.         jsr     fmpy_xy                 ;multiply x*y first
  82.         move    a,x1                    ;save product mantissa
  83.         move    b,x0                    ;save product exponent
  84.         move    fp_space:fp_temp+1,a    ;restore ma
  85.         move    fp_space:fp_temp+2,b    ;restore ea
  86.         jmp     fsub_xa                 ;subtract product from a
  87.  
  88.